Bloch oscillations in complex crystals with VT symmetry 
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Bloch oscillations (BO) in complex lattices with VT symmetry are theoretically investigated with 
specific reference to optical BO in photonic lattices with gain/loss regions. Novel dynamical phe- 
nomena with no counterpart in ordinary lattices, such as non-reciprocal BO related to violation of 
the Friedel's law of Bragg scattering in complex potentials, are highlighted. 
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Bloch oscillations (BO), i.e. the coherent oscillatory 
motion of a quantum particle in a periodic potential 
driven by an external dc force, represent one of the most 
striking predictions of wave mechanics of periodic sys- 
tems. In the quantum realm, BO have been experimen- 
tally observed for electrons in semiconductor superlat- 
tices [H and for ultracold atoms or Bose-Einstein con- 
densates in tilted optical lattices [2j] . Classical analogues 
of BO have been also proposed and observed for optical 
3 and acoustical ^3] waves. The phenomenon of BO is 
related to the transition of energy spectrum from contin- 
uous (energy bands) to (nearly) discrete, with the forma- 
tion of Wannier Stark (WS) ladders when the dc force is 
applied. In physical space, BO are generally explained 
in terms of wave Bragg scattering off the periodic poten- 
tial, causing a wave packet to oscillate rather than trans- 
late through the lattice. Previous studies on BO have 
considered particle motion in real- valued potentials, may 
be in presence of nonlinearity, lattice disorder or parti- 
cle interactions [1]. However, particle or wave dynamics 
in complex potentials is not uncommon in physics 
1^. Complex potentials for matter waves emerge in the 
near resonant interaction of light with open two-level sys- 
tems 0, whereas in optics wave propagation in lossy or 
active waveguides is described by a complex refractive 
index Q. Complex potentials have received in recent 
years an enormous interest since the discovery that non- 
Hermitian Hamiltonians with parity-time {VT) symme- 
try may possess an entirely real-valued energy spectrum 
below a phase transition (symmetry-breaking) point jTo| . 
VT Hamiltonians belong to the more general class of 
pseudo-Hermitian systems [TT'|. Optical realizations of 
T'T-symmetric potentials, based on light propagation in 
coupled waveguides with a complex refractive index, have 
been recently proposed p^ - [T^ , whereas energy band for- 
mation and spontaneous symmetry-breaking thresholds 
have been investigated in a few VT potential models [l^ 
and in presence of nonlinearity [16[ or lattice disorder 
It is the aim of this Letter to investigate the onset 
of BO in periodic lattices with VT symmetry, highlight- 
ing some unique features arising from the complex nature 
of the scattering potential. We consider here BO for light 
waves in photonic lattices 's'l, however the analysis can 
be applied to other physical systems, such as to matter 
waves in accelerating complex optical lattices. 



Paraxial propagation of light at wavelength A in a 
weakling-guiding one-dimensional waveguide lattice with 
a superimposed transverse refractive index gradient is de- 
scribed by the Schrodinger-type equation [11, [3| 



iXd^ij ^ -—dl^p + V{x)ip + Fxip = {-Ho + Fx)^ (1) 
2ns 

where A — X/{2tt) is the reduced wavelength, ris is the 
substrate refractive index, V{x) = Ug — n(x) is the po- 
tential, 71(2;) = n(x -|- a) is the refractive index profile 
of the lattice (spatial period a), and F is the refractive 
index gradient which mimics an external dc force. For 
F — VT symmetry requires V{—x) = V*{x), i.e. suit- 
able combinations of optical gain and loss regions in the 
lattice, as discussed in |12h14| . The real and imaginary 
parts of the potential are denoted by Vr{x) and aVi{x), 
respectively, where a > is a dimensionless parameter 
that measures the anti-Hermitian strength of Hq. For 
given profiles of Vr and V7, the spectrum of 'Hq turns 
out to be real- valued and composed by bands separated 
by gaps for a < ac, where ac ^ corresponds to the 
spontaneous symmetry breaking point. For a > ac, band 
merging and appearance of pair s of complex-conjugates 
eigenvalues is observed [Tsl Il5j. For instance, for the 
potential 

Vr{x) = yocos(27ra;/a), Vi{x) = sin(27ra;/a) (2) 

one has ac — I (Tsj . 

Let us first consider the case a < ac, and indi- 
cate by 4>n{x,K) — Un{x, k) expiiKx) the Bloch-Floquet 
eigenfunction of Ho, corresponding to the n-th lattice 
band and to wave number k in the first Brillouin zone 
(— 7r/a < K < 7r/a) with periodic part Un{x,K). Thus 
Hoclinix, k) = En{K)(j){x, k), whcrc £"„(«;) is the real- 
valued band dispersion curve. As for real potentials, 
one can show that £"„(—«) = £"„(«;). Additionally, 
Re(u„(K,x)) and Im(u„(K,a;)) have well defined and op- 
posite parity under the inversion x —x, whereas 
(j>l-^{x,K) = 4>n{—x,~Hi) are eigenfunctions of the ad- 
joint "HJ with same eigenvalues £"„(«). Assuming that 
/„ = dxu'^{—x, —K)un(x, k) does not vanish, i.e. in 
absence of spectral singularities, the orthonormal condi- 
tions {(t>l{x,K')\(j)n{x,K)) = J^^dx(l)l*{x,K')(j)n{x,K) = 

V>nSn.iSii^ — 1^') hold, where ~ ±1 is given by the 



sign of /„ in each band [Ifl, ■ To investigate WS lo- 
calization, like in ordinary crystals we expand the wave 
packet ■0 as a superposition of Bloch-Floquet states [l^ , 
z) = J2n I dKCn{K, z)(j)n{K, z). The scalar product of 
both sides of Eq.(l) with (0]^(a;, k)| yields the following 
evolution equations for the spectral coefficients c„(k, z) 



where Xn^i{n) = {2TTi/a) dxu'^{—x,—K)dK.ui{x,K). 
The off-diagonal elements Xn,i [n ^ I) in Eq.(3) 
are responsible for interband transitions, i.e. 
Zener tunneling (ZT), and can be written in 
the equivalent form Xn^i — (2ttX'^ /ans)[En{K) — 
Ei{k)]~'^ u'^{—x, —K)dxV{x)ui{x, k). If bands n 
and / are separate by a large gap and the force F is 
small enough such that F\Xn,i{K)\ <^ \En{K) — Ei{k)\ 
in the entire Brillouin zone , ZT is negligible as in 
ordinary lattices [l^. Taking Xn,i ~ for n 7^ Z 
(single- band approximation), the WS spectrum of 
the n-th band is then obtained by letting in Eq.(3) 
Cn{K, z) = a„(K) exp(— i£'z/A) and imposing the 
boundary condition a„(27r/a) = a„(0). This yields 
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FZa +(£;„)+ ^7„ G-0,±l,±2,.. 



where we have set {En) = (a/27r) /_ 
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(4) 



and 



7„ = (2?„i^a/27r) J^^J^,^dKlm{Xn,n{K))- As compared to 
an ordinary crystal [19|, the main difference here is that 
the WS ladder spectrum is complex- valued, even though 
a < ttc- This is simply due to the fact that the external 
force F breaks the VT symmetry of the Hamiltonian. In 
each ladder, the imaginary parts of the eigenvalues have 
the same value (i.e., independent of /), equal to 7„. Since 
Im(X„^„(K)) is an even function of k, 7„ is usually non- 
vanishing, and its sign is reversed as the direction of the 
force, i.e. the sign of F, changes. Therefore, in complex 
lattices below the phase transition point, BO with spa- 
tial period zb = X/{Fa) are observed, however they are 
either amplified or damped depending on the sign of F. 
This behavior is illustrated in Fig.l, which shows the on- 
set of amplified or damped BO as obtained by numerical 
simulations of Eq.(l) for the complex lattice of Eq.(2) 
with a — 0.3 and for an input excitation of the lattice 
with a broad Gaussian beam at normal incidence. In the 
figures, the evolution of the total beam power P — P{z) 
is also shown. Note that, in spite all WS eigenvalues have 
the same imaginary part, P{z) shows a non-exponential 
behavior, with abrupt changes of P{z) in correspondence 
of Bragg scattering of the beam from the lattice. This 
behavior is ascribable to the non-orthogonality of WS 
states in complex potentials. 

As the non-Hermitian parameter a is increased to 
cross the phase transition value Oc, gap narrowing till 
band merging, associated to the appearance of pairs of 
complex-conjugate eigenvalues, is observed 13[, making 
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FIG. 1: (color online) BO in the complex lattice of Eq.(2) for 
a = 0.3 and for an index gradient (a) F — —3.52 m^^, and 
(b) F — 3.52 m~^. Other parameter values are a — Q fim, 
A = 633 nm, Vb = 0.001, and = 1.42. The figures show 
the evolution of field amplitude |'i/'(a;,z)| and beam power 
P = / dx\4>{x, z)\^ (normaUzed to the input power) versus 
propagation distance z. The BO period is = 3 cm. 



the previous analysis inapplicable. Here we will limit to 
investigate the onset of BO in the a = ac and a > ac 
regimes by considering two specific examples which en- 
able an analytical treatment. 

The first example is given by the potential (2) at 
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FIG. 2: (color onhne) Same as Fig.l, but for a = Oc = 1 and 
(a) F = -7.03 m~\ (b) _F = 7.03 m~\ Other parameter 
values are as in Fig.l, except for Vb = 2 x 10~*. Note that 
BO disappear for a negative index gradient F. 

the phase transition point a — ac = 1, i.e. V{x) — 
Vq exp{2Trix / a) [l^], for which band gaps disappear and 
the energy dispersion is described by the free-particle 
parabola E{k) — (A^/2ns)K^, periodically folded inside 
the first Brillouin zone. In presence of an external force 
F, numerical simulations indicate that the wave packet 
evolution strongly depends on the sign of F, as shown 
in Fig. 2. Remarkably, BO disappear for F < 0, and the 
wave packet propagates following a parabolic path as if 
the complex lattice were absent [Fig. 2(a)]. This result in- 
dicates that a net drift along the lattice can be achieved 
by the application of a dc force, an effect which is not pos- 
sible in ordinary crystals without collisions. Conversely, 
for i<" > an oscillatory dynamics survives, with spatial 
period equal to zb = X/{Fa) [Fig. 2(b)]. In this case, 
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at propagation distances corresponding to Bragg reflec- 
tion from tlie lattice, a rather abrupt amplification of 
the scattered beam is observed [see the inset of Fig. 2(b)]. 
This dynamical regime, which can be referred to as uni- 
directional BO, has no counterpart in ordinary lattices 
and is related to violation of Friedel's law of Bragg scat- 
tering by complex potentials 0, H^l, i-e. breakdown of 
Bragg diffraction invariance by crystal inversion. Ana- 
lytical insights into unidirectional BO can be gained by 
writing Eq.(l) in momentum space and assuming that 
the lattice is excited by a broad wave packet at normal 
incidence, i.e. ip{x,0) — J dkg{k) exp{—ikx), where the 
plane-wave spectrum g{k) is narrow at around k = 
with a spectral width much smaller than the Bragg wave 
number fcs = lix j a. Let us indicate by ipoix, z) the prop- 
agated wave packet that one would observe in absence of 
the lattice, i.e. the solution to Eq.(l) with y = 0. The 
center of mass {x){z) of ?/'o(a^:'Z) follows the parabolic 
path {x){z) — — (F/2ns)z^ due to the external force F, 
whereas the mean value of momentum {k) grows lin- 
early, {k){z) = —Fz. In presence of the complex lattice, 
Bragg diffraction sets in when (k) gets close to —ks/i, 
but not when (k) ~ fcB/2; this nonreciprocal behavior 
of Bragg scattering by complex potentials is discussed 
e.g. in p, [23|. In presence of the lattice, one can show 
that the wave packet is expressed by a superposition of 
retarded wave packets ipo{x,z — nzs) with amplitudes 
Bn{z), namely 



ip{x,z) 



+ 00 

E 

n=0 



B„{z)iPq{x, z - hzb). 



(5) 



In Eq.(5), Bq{z) = 1, whereas for n > 1 the amplitudes 
Bn can be iteratively calculated from the relation 



Bn{z) 



■1^ 



d^exp[-iipn{0]Bn-iiO, (6) 



where we have set (/j„(z) {X'^kB/2nsF)[{Fz/X - 
kBu)"^ + ksiFz/X - ksn) + fc|/3] . In the limit zb > 
{nsa/irFy^^ and for F < 0, \Bn>i{z)\ remains small be- 
cause of ipn{z) has no stationary points for z > 0. Physi- 
cally, this corresponds to the absence of Bragg scattering 
(the condition (fc) ^ —kB/2 is never satisfied), and the 
scenario of Fig. 2(a) is observed. Conversely, for F > 
the phase has a stationary point at z = ZB{n — 1/2), 
and one obtains |i?„(z)| ^ r"i?(z — uzb + zb/2), where 
r = (Vb/A)(nsa/F)i/2 and H{^) is the Heaviside step 
function [H(0 = for ^ < 0, H{£_) = 1 for ^ > 0]. 
In this case, Eq.(5) clearly shows that a periodic regen- 
eration and amplification of the wave packet ipo due to 
Bragg scattering occurs at z = zb /'2,3zb/'2,5zb/2, 
explaining the behavior observed in Fig. 2(b). For param- 
eter values used in the simulations, the predicted value 
of amplification factor F is F ~ 2.186, with which the 
staircase behavior of beam power shown in the inset of 
Fig. 3(b) is reproduced with excellent accuracy. 

As a second example, we consider a complex binary 
lattice [l6j with broken VT symmetry at ac = 0, and 



show that the external force can restore a real-valued 
energy spectrum. The lattice consists of an array of 
equally-spaced single-mode waveguides with alternating 
gain (V7 > 0) and loss {Vj < 0) channels, i.e. Vr{x) — 
Y,i h{x -la- a/2) and Vi{x) = Y.i{-^)^Kx " - a/2), 
where h{x) = h{—x) defines the (real) index profile of 
the waveguide channel (see Fig. 3(a)). Note that, for 
a ^ 0, the lattice period is 2a. A typical band diagram 
is shown in Fig. 3(b). In the tight-binding approximation 
and for a ^ 1, light transport in the lattice is described 
by coupled-mode equations for the modal aniplitudes ai 
of light trapped in the waveguides 0, [H, [H, [2l| 

" ~t("'+i + '''^ + + ^(-l)'f «/ (7) 
where A/4 is the real- valued coupling rate between ad- 
(a) .(b) 
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FIG. 3: (color online) BO in a complex binary lattice (a = 
10 /im). (a) Profiles of Vr and Vi, and (b) corresponding band 
diagram for a = 0.01 (real and imaginary parts of energy E). 
(c) Resonance curves in the (A/ f, g/ f) plane corresponding 
to a real-valued WS ladder spectrum, (d) Breathing BO mode 
(single waveguide input excitation) for F = 1.618 m^^, cor- 
responding to a Bloch period zb — 39.1 mm. 

jacent waveguides, / — Fa/X, and g is the real- valued 
power gain/loss coefficient in the waveguides, which is 
proportional to a. For the potential of Fig. 3(a) with 
a = 0.01, from the band diagram of Fig. 3(b) (the two 
merged bands n — 1 and n = 2) one can estimate 
g ~ 1.46 cm^^ and A ~ 14 cm"^. For F = 0, the 
spectrum of Eq.(7) is purely continuous and composed by 
two merged minibands [the two lowest bands of Fig. 3(b)] , 
namely [16| 



E{k) = ±Av/(A/2)2 cos2(Ka) - {g/2)^ , |k| < n/2a. 

(8) 

Since complex-conjugate eigenvalues appear near the 
band edges for any value of 5 > 0, the lattice has zero 
threshold for VT symmetry-breaking, i.e. Uc = 0. In 
presence of the external force F, the spectrum becomes 
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purely discrete and composed by two interleaved WS lad- 
ders eI^'' and as in a binary lattice with real po- 
tential [2l|. Specifically, one can show that 

i^a.^) ^ 2Fal - + (9) 

TT TT 

= 0, ±1,±2,...) where yi,2 = |?/i,2| exp(i.^i^2) = 
Re(cr) ± [1 -HRe(cr)2]i/2^ ^ = ^^^^^ g^^^^j the 2 x 2 matrix S 
is the propagator for the differential system dxi^2/dK ~ 
±i(A/4/) cos(K/2)a;i,2 -I- ig/4:f )x2A from k = to k 
27r. Note that the two WS ladders are generally complex- 
valued with opposite sign of their imaginary part. How- 
ever, if the condition Re{a) = is satisfied, yi,2 = ±1 
and the spectrum becomes real-valued and given by a 
single WS ladder, namely Ei = Fal {I = 0, ±1, ±2, ...). 
Figure 3(c) depicts the numerically-computed locus in 
the (A//, g//) plane where Re{a) = 0, showing the exis- 
tence of a sequence of resonance curves. A perturbative 
calculation of the propagator S in the g — ^ limit shows 
that the resonance curves depart from the points of the 
A// axis that are the roots of Bessel function Jg [2ll |. 
Along such resonance curves, a forced Bloch particle un- 
dergoes undamped nor amplified BO, with spatial period 
zb — X/{Fa). This is shown in Fig. 3(d), which depicts 
an example of a BO breathing mode, obtained by nu- 
merical analysis of Eq.(l) for single waveguide excitation 
of the complex binary lattice of Fig. 3(a), with a forcing 



F = 1.618 m~^, corresponding to point A {g/f ~ 0.91, 
A// ~ 8.72) on the third resonance curve of Fig. 3(c). 
Photonic VT lattices may provide an accessible labo- 
ratory to experimentally observe BO in complex crys- 
tals. Complex refractive indices could be realized in en- 
gineered arrays of active semiconductor waveguides with 
selective pumping to achieve gain and loss regions. At 
optical or near-infrared wavelengths, a modulation depth 
of the imaginary refractive index An/ ^ 2 x 10^** (as in 
Fig. 2) requires a gain/absorption coefficient of the order 
^ 10-20 cm~^, which is available with current semicon- 
ductor optical amplifier technology. Smaller gain/loss 
coefficients are needed for the example of Fig. 3. Non- 
reciprocal BO could be also observed in passive (i.e. with- 
out gain) waveguide arrays with loss modulation, in a 
geometrical setting similar to the one recently realized 
to demonstrate VT symmetry breaking in a passive Al- 
GaAs optical coupler [22l | . The transverse index gradient 
can be realized by either chirpin g th e waveguide width 
Q or by waveguide axis bending |l8| . 
In conclusion, BO in complex crystals with VT symmetry 
behave differently than in ordinary lattices owing to non- 
reciprocity of Bragg scattering. Our results indicate that 
the conventional wisdom of BO and transport for classi- 
cal or matter waves needs to be modified when dealing 
with complex crystals. 
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